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Majorana Edge Modes of Superfluid ^He A-phase in a Slab 
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Motivated by a recent experiment on the superfluid ^He A-phase with a chiral p-wave pairing 
confined in a thin slab, we propose designing a concrete experimental setup for observing the 
Majorana edge modes that appear around the circumference edge region. We solve the quasi- 
classical Eilenberger equation, which is quantitatively reliable, to evaluate several observables. 
To derive the property inherent to the Majorana edge state, the full quantum mechanical 
Bogoliubov-de Gennes equation is solved in this setting. On the basis of the results obtained, 
we perform decisive experiments to check the Majorana nature. 
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There has been growing interest in Majorana quasi- 
particles (QPs) and fermions, which are expected to play 
a major role in various research fields, ranging from cos- 
mology and high-energy physics to condensed matter 
physics,^-' and even to topological quantum computing.^-' 
The Majorana QP and fermionic operator are defined 
by 7 = 7^ and 'J = ^'t, respectively, which imply that 
the particle and antiparticle are identical and thus are 
electrically neutral. They are an intriguing subject to 
further study in their own right. It has been proposed 
that the Majorana nature brings new physics, such as 
the non-Abelian statistics for spinless chiral superfluids^' 
and Ising-like spins for time-reversal invariant superflu- 
ids.^"^-* Obviously, we like to have more candidate mate- 
rials to realize it. 

The zero-energy bound state appears whenever the un- 
derlying potential for QP changes its sign. This is the so- 
called TT-phase shift physics exemplified by charge or spin 
density waves, ^'''^ FFLO states,^' or stripes in high-Tc 
cuprates."'^'^'' The Majorana QP is a special case of this 
type: The Majorana conditions'^) under which the Majo- 
rana QP 7q = 7o is produced necessarily at a zero energy 
are summarized as (1) chiral p-wave superconductors or 
superfluids where the Bogoliubov QP 7]^ with an energy 
E satisfies — ^~^e owing to particle-hole symmetry, 
and (2) a certain type of vortices with an odd winding 
number. For the edge or surface, the Majorana fermion 
\E' = \E'1' exists as well as =70.^'"'^' 

There is as yet no firmly established chiral p-wave su- 
perconductors so far, although several theories'''''^''®-' 
have explored its possibility in Sr2Ru04. We have to re- 
sort to the use of the superfluid '^He system in which we 
may manage to manipulate the order parameter (OP) 
so that the Majorana conditions are fulfilled by setting 
up an appropriate boundary condition. A good exam- 
ple of this kind has been recently presented by several 
groups. ■^"^^ They demonstrate that the Majorana edge 
modes in "^He-B give rise to the anisotropic spin dynam- 
ics, while the OP in the B-phase is a time-reversal 
invariant, which is in contrast to that in the A-phase. 
Thus the main purpose of the present work is to propose 



a concrete experimental design for detecting the Majo- 
rana QP based on quantitatively accurate microscopic 
calculations for the '^He A-phase. 

Here, we examine the superfluid ^He A-phase, which is 
characterized by a chiral p-wave type OP without doubt 
which is a time- reversal breaking state. '^•' Thus, the Ma- 
jorana QPs in the A- and B-phases are worthy of com- 
parative studies. In fact, as we will see later, we can treat 
the Majorana QPs in the A- and B-phases in the same 
experimental setup in a slab geometry by merely chang- 
ing temperature or pressure, which is a large advantage 
over other proposed systems. The system provides a use- 
ful testing platform for checking various ideas associated 
with the so-called topological order. 

Our basic idea is stimulated by the experimental setup 
used by Bennett et al.,^^'^ where a superfluid '^He A- 
phase is conflned in a thin slab with a thickness Z? = 0.6 
/im and {L — 10mm) x 7mm, as schematically shown 
in Fig. 1. Other experimental groups "'^^"^'^ use various- 
thickness samples, ranging from D = Q.l to 10 /im. Since 
the so-called dipole coherence length ^rf ^ 6 /im,'^' the 
/-vector, which characterizes the chiral direction of the 
A-phase OP, points perpendicular to the slab's upper 
and lower surfaces (B), as shown in Fig. 1. We assume 
the specular boundary condition in this paper. The four 
sides of the rectangular strips of the slab are an ideal 
place to accommodate the Majorana edge modes. 

We use two theoretical methods, the quasi-classical 
Eilenberger equation and Bogoliubov-de Gennes (BdG) 
equation. The quasi-classical framework is used for de- 
tailed studies of the quantitative aspects of the super- 
fluid properties at the edges. Using the BdG framework, 
we reveal the Majorana nature and quantum-mechanical 
spectral structures of the low-lying edge modes. 

We consider the following spatially one-dimensional 
problem along the a;-axis with a length L [0 < x < L). 
See section A in Fig. 1 where, along the z-direction, the 
system is uniform because the i-vector points perpen- 
dicular to the upper (B) and lower faces. The length is 
scaled by the coherence length at T = 0, ^0 ~ 80 nm. In 
this situation, the pair potential of the chiral OP in the 
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Fig. 1. (Color online) Schematic configuration of the whole sys- 
tem. A: the cross-section considered here. B: Upper surface of the 
container. The Z-vector is locked to the upper and lower surfaces. 
The mass current j flows along the edges. 

A-phase is described by 

A{x,k) = rix{x)kx: + r]y{x)ky. (1) 

Note that we are considering the three-dimensional Fermi 
sphere for '^He in the momentum space. The relative 
momentum of the Cooper pair on the spherical Fermi 
surface is given by fc = (sin 6k cos (jjk , sin 9k sin (j)k , cos 9k ) , 
and the Fermi velocity is given hy v = vpok. It is im- 
portant to notice that the quasi-particles described by 
Green's functions move under the "given" pair potential 
A(x, fc). The relative phase between rjx and r]y is 7r/2, 
since the i-vector points to the 2;-direction. 

The applied field H is an important controlling param- 
eter of the system: For the slab thickness -D < ^d, the field 
perpendicular (parallel) to the slab should he H>Ha {H 
arbitrary) so as to satisfy the condition that dLH with 
being a dipolar field ^ 2mT.^^^ Then, an equal spin 
pairing is realized and we can always ignore the + 
Cooper pair by changing the spin quantization axis ap- 
propriately.^'^) Thus, we omit the spin degrees of freedom 
in our calculation, e.g., in eq. (1). 

It is known that different from that in the bulk "^He, 
in the slab of superfluid ■^He, the A-phase is stabilized 
at much lower temperatures T down to the pressure 
P = 0.24) rpj^jg condition is fulfilled for 0.06 < £> < 0.3 
/im where the A-phase is stable at T = and P — {). For 
sub-/im D the A-phase changes into the B-phase with de- 
creasing T . Thus, it allows us to treat both phases under 
the same experimental setup. A-B control is also possi- 
ble by varying pressure. Theoretically, our calculations 
are reliable quantitatively because superfluid '^He in the 
lower-pressure regions of interest here can be described in 
a weak coupling scheme without delicate strong coupling 
effects.25) 

We start with the quasi-classical Eilenberger equa- 
tion, ■^^"•^'') which has been used in the study of '^He 
superfluidity. '^^"■^'^) The quasiclassical Green's functions 
g{uJn,k,x), /(w„,fc,x), and f''{u!n,k,x) are calculated 
using the Eilenberger equation 

{u,, + v-V)f^Ag, {u„~v-V)P =A*g, (2) 

where (1 — ff^y^'^, Keg > 0, and v — v/v-po- We solve 
eq. (2) by the Riccati method, '^4''^^'' with the specular 
boundary condition at both edges: a; = and L. The 
temperature is measured using the transition tempera- 
ture Tc. Matsubara frequency a;„ = {2n+ l)7rfcBT, energy 
E, and pair potential A are in a unit of irk^Tc. 
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Fig. 2. (Color online) Order parameter profiles along the x-axis 
for (a) L = 100 and (b) L = 12 in a unit §o. 

The order parameter rij{x) (j — x,y) is self- 
consistently calculated by 

i^,{x)^ig^N^T (^^^(^ + -^^*))fc' 

0<cj„<ti;cut 

with {goNo)-^ = lnr + 2TEo<.„<..„,Kr^ {■■■)!. 
indicates the Fermi surface average, and iVo is the density 
of states (DOS) at the Fermi energy in the normal state. 
We use Wcut = 40A:b?c. The temperature is flxed at T = 
0.2Tc throughout the paper. 

Using the obtained self-consistent solutions, the mass 
current is calculated using 

ji^) ^ ijx,jy,jz) o^T {vling)k. (4) 

When we calculate the quasiparticle states, we solve 
Eq. (2) with ia;„ E + i5. We typically use S — 0.005. 
The local density of states (LDOS) for quasiparticles nor- 
malized by iVo is obtained as 

N{E,x) = {K,e{g{u!n,k,x)\i^^^E+is})k- (5) 

In Fig. 2(a), we show the OP profiles for L= 100. It 
is seen that, because of the specular boundary condition 
at the edges, the component, which is perpendicu- 
lar to the edges, becomes zero, while the parallel com- 
ponent ky is enhanced by compensating for the loss of 
the kx component. Thus, at the edges, the polar state 
is realized. Towards the center, the kx {ky) component 
increases (decreases). The A-phase with kx + iky is at- 
tained in the central region where the ^-vector points to 
the z-direction. We also show the result for L = 12 in 
Fig. 2(b), where even in the central region the complete 
A-phase OP is not recovered and the polar phase nature 
is dominant there. 

Figure 3(a) shows LDOS at the edge x — and the 
center x — bQ. It is clearly seen from curve 1 for a; = that 
there exists a zero-energy state with a substantial weight 
at i? = corresponding to the Majorana edge mode, that 
is, LDOS is expressed as N{E, x = 0) = 70 + aE^ near 
E — Q. The first (second) term comes from the Majorana 
edge mode (point node contribution of the chiral state 
in the bulk A-phase). LDOS at a; = 50 (curve 2) shows a 
typical point node behavior of N{E,x = 50)(xE^. 

In Fig. 3(b), we show the extension of the Majo- 
rana edge mode towards the center from the edge at 
a; = 0, which spreads out an order of ~5^o- In Fig- 3(e) 
the spectrum evolution from the edge to the center is 
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Fig. 3. (Color online) Calculated results for L = 100. (a) LDOS 
N(E,x) at (1) x = and (2) x = 50. (b) Zero-energy LDOS 
N{E = 0, x) as a function of x from the edge at a; = 0. (c) Spectral 
evolution N{E,x). The cross sections at x = and x = 50 are 
same as in (a) and _E = as in (b). (d) The angle-resolved LDOS 
N{E, X, Sfc) in the momentum space at a; = with <pk = and 
6»fe = 89.5° (1), 60.5° (2), 45.5° (3), 30.5° (4), 15.5° (5), and 0.5° 
(6). The zero-energy states are accumulated from all flj, values, 
except for (6). (e) Mass current jy{x). 

shown. The spectrum of the Majorana edge mode grad- 
ually changes into the bulk spectrum. The angle-resolved 
LDOS N{E, X, Ok) is shown in Fig. 3(d), where Ok is the 
polar angle from the fc^-axis of the 3D Fermi sphere. 
The zero-energy LDOS comes from the quasi-particles 
with the nonvanishing component, meaning that the 
QPs reflected by the edge form the Andreev bound state 
exactly at i? = because the QPs effectively feel the 
sign-reversed pair potential upon retracing the incoming 
path. The physics of the 7r-phase shift is working here. 
The mass current jy{x) is displayed in Fig. 3(e), show- 
ing that the chiral current flows circularly along the four 
edges, as shown in Fig. 1. By applying a magnetic field 
either parallel or perpendicular to the slab, we can pro- 
duce a spin imbalance due to the Zeeman shift between 
the 'l^t pairs and the ]^\. pairs. This results in a net spin 
current in addition to the above mass current. 

The LDOS's at the edge are calculated for various 
lengths in Fig. 4(a). As L decreases, the interference be- 
tween two edges, which makes the zero energy levels split, 
becomes important and results in a decrease in the spec- 
tral weight of the zero-energy Majorana QP. The spectral 
weight of the Majorana QP N{E = 0, a; = 0) is shown as 
a function of L in Fig. 4(b). Beyond L ~ 20, its weight 
remains constant. 

To examine the full quantum nature of the discretized 
Majorana levels, we solve the BdG equation'^^"'^*' with 
the pair potential obtained by the Eilenberger equation 
(2). Let us consider the situation with d\\ z for exam- 
ple. The BdG equation can be separated into two spin 



(a) (b) 




Fig. 4. (Color online) (a) LDOS at the edge a; = for (1) L = 100, 
(2) L = 12, and (3) L = 8. (b) The spectral weight of the zero- 
energy state N{E = 0,x = Qi) as a function of L. 

sectors labeled a — ±, which denotes the eigenstates of 
the Pauli matrix fj,, corresponding to the choice of the 
spin quantization axis parallel to y. Then the BdG equa- 
tion reduces to each spin sector of the QP with the wave 
function ipn.a = [un.a,Vn,aV and the energy Eny.^^^^^'> 

j dr2'H[ri,r2)<Pn,<y[r2) = En(Pn.a{ri). (6) 

Here, we set = - {1-0)22 = - ^ - and (7i)i2 = 

^('"i, '"2) = ~('^)2i- The eigenstates of eq. (6) yield one- 
to-one mapping between the positive energy states ipE 
and the negative energy states y-_E — owing to the 
symmetry % — —Tx'H*fx, leading to the relation of the 
Bogoliubov QP operator 7b =7! ^j, as mentioned above. 

Before turning to the numerical results of eq. (6), 
we should mention the Majorana nature of the edge 
states. Within the weak coupling regime fcp^^l, eq. (6) 
with eq. (2) is solved for the edge states as ipk,a = 
e'^*'/'^e'^-^\if{k^,x) with the energy E^-Aq^, where 

K={ky,k,^), /(fc,x)=sin(fca;)e-^/««, and i9 is the U(l) 
gauge of A. Then, the field operator in the y-quantization 
axis is expanded in terms of the Bogoliubov operators 
IK as = Efc[e"-^""7i^ + h.c.]/(fc,,a;)e'''/2 + ^-tuik, 
where ^'buik denotes the contribution from the bulk ex- 
citations. Although 5'buik contains the low-energy states 
due to point nodes, their contributions are negligible if 
T <C Tc is considered. This is because the DOS due to the 
edge states yields N{E) ocE'" near E^E-p, which over- 
whelms the DOS due to the point nodes N{E)(xE'^, as 
shown in Figs. 3(a) and 3(b). This predicts the Majorana 
condition, e"-^ '^l. for T and 7if = 7!^^. 

In Figs. 5(a) and 5(b), we show the spectrum of 
the Majorana edge modes En{K) normalized by Aq at 
r = 0. Here, the nonlocal pair potential in eq. (6) is ob- 
tained by the self-consistent solution of the Eilenberger 
equation (2) through A(ri,r2) = / ^^A(r, fc)e^'^ ''i2, 
where ri2~ri — r2 is the relative coordinate. Along the 
fca;-axis, En{K — 0) is dispersionless and along the ky- 
axis two ky-lmeai dispersions near E — appear, each 
coming from the left and right edges. At ky — 0, the two 
Majorana modes are situated exactly at E — 0. A stere- 
ographic view of those modes is shown in Fig. 5(c). It is 
seen that the linear Dirac dispersion continues along the 
ky direction, forming a "Dirac valley" . The Dirac valley 
merges into point nodes at k^ — ikp- These results are 
consistent with those obtained using the previous Eilen- 
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Fig. 5. (Color online) Results of BdG calculations at T = 0.2Tc, 
where A{x,k) is taken from Fig. 2(a) with L = 100. Eigenvalues 
of the low-lying excitations along (a) kz for ky =0 and (b) ky 
for kz = 0. "R"("L") for right (left) edge, (c) 3D view of the 
dispersion relation: "Dirac valley" . 
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tice of electrons placed on the surface. We need a special, 
but feasible configuration of the experimental setups. 

In summary, we have designed a concrete experimen- 
tal setup to observe the Majorana particles at the edge 
in a certain slab geometry. We calculate the microscopic 
Eilenberger equation to yield quantitatively reliable in- 
formation on observable physical quantities in realis- 
tic situations for the ^He A-phase. The Bogoliubov-de 
Gennes equation is solved to explicitly demonstrate the 
Majorana nature of edge states. Several feasible and ver- 
ifiable experiments for checking the Majorana nature are 
proposed. 

The authors thank K. Kono and J. Saunders for infor- 
mative discussions on their experiments. 
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berger equation and reveal the detailed quantum level 
structures of Majorana modes. 

There are several experimental ways to detect Majo- 
rana modes. The most direct evidence of the Majorana 
nature is the observation of the anisotropy of spin sus- 
ceptibility. Using the Majorana nature of the edge states 
vl/j_ = e*''\l/j_ for T < Tc, as derived above, the local 
spin operators result in Sx ~ Sy ^ and only the z- 
component parallel to the d-vector remains nontrivial, 
5*2 « — 2i4'+\['_. This predicts the Ising-like spin dynam- 
ics in the A-phase as well as in the B-phase.'^'^-' This is 
in sharp contrast to the anisotropy in the bulk A-phase, 
where the susceptibility parallel to d is suppressed at low 
T according to the Yosida function. 

The surface specific heat measurement, which was per- 
formed previously^^^ in connection with the Andreev sur- 
face bound state detection in the B-phase, resolves its 
contribution Csurface(T) =7sT because 7s ociV(i? = 0, x = 
0) from the bulk contribution with Cbuik cx coming 
from the point nodes where N{E) oc E^. Note that, in 
the slab geometry, there is no low-lying excitations other 
than C oc T'^ from the two upper and lower surfaces (B 
in Fig. 1). Thus, the Csurfaco(T) = jgT contribution of 
the Majorana QPs is distinctive. The interference effect 
between Majorana QPs when L decreases, each coming 
from both edges, yields the LDOS N{E)^N{0)+a\E\, as 
shown in Fig. 4. This extra linear \E\ term gives the spe- 
cific heat CsurfaccCT") (xT^. The relative weight of these 
T-linear and terms depends on the distance L be- 
tween the two edges, which is precisely evaluated (see 
for example L = 12 case (2) in Fig. 4(a)). 

Quasi-particle scattering or QP beam experiments are 
extremely interesting. They were performed in the past 
on ^He where roton-roton scattering is treated^°^ and on 
the ^He B-phase where the Andreev surface bound state 
is investigated. ^^-^ Using this method, we may pick up 
Majorana QPs with a particular wave number because 
we like to manipulate Majorana QPs located at ky = 0, 
which is separated from other QPs in the nodal region 
at kz = ±kp. 

The other option might be to use a free surface where 
the Majorana surface bound state is formed. As shown 
by Kono's group, ''^^ the bound state can be detected 
through the excitation modes of the fioating Wigner lat- 
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